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ABSTRACT 


This investigation concerns with the problem of 
a shell of translation, shallow in the longitudinal 
direction, acted upon by uniformly distributed vertical 
load and wind loads. A unit of a shell is formed by * 
translating an inclined line vertically over a parabola. 
Sucfr units can be combined to form an ’arched folded 
plate*, the cross section of which is of a conventional 
folded plate where as the longitudinal axis is a parabolic 
cui've. The transverse edges are either fixed or hinged. 

The longitudinal edges have toeen kept free, hinged or 
fixed. 


The equilibrium equations, stress -strain laws and 
strain-displacement relations are derived according to 
linear theory of shells. The equilibrium equations are 
reduced to three equations in terms of fundamental displaoe 
mants. Assuming the nature of displacements in the longitu- 
dinal direction, the unknown functional values in the trans 
verse direction are evaluated by Runge-Kutta Gill integra- 
tion procedure which yields the stresses and displacements 
in the structure. Several numerical results are presented. 
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CHAPTER 1 


INTRODUCTION 

The economy of shell or folded plate structures 
is based on the use of main structural system as the 
building enclosure* A wide range of shells is made use 
of in Aeronautical, Civil and Naval structures which 
necessitates a high degree of sophistication in analysis., 
Most of the shells used for roofing comes under the 
category of thin shells where the ratio of deflection 
to thickness remains much smaller than unity. 

The method of analysis of shells and plates 

under the 'thin shell theory’ can be carried out on the 

* 

assumptions proposed by Kirchoff |l|; viz. (a) the 
straight fibres of the plates perpendicular to the 
middle surface before deformation remain so after defor- 
mation and do not change their length, (b) the normal 
stress acting on planes parallel to the middle surface 
can be neglected in comparison with the other stresses 
without loss of much accuracy. This was elucidated by 
Novozhilov |2| who proved that the errors introduced by 
Kirchoff s assumption are of the order of h/R compared 
with unity. The relation between forces, moments and 


* Numbers within | | designate references at end of thesis. 
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displacements at the middle surface are accepted with 
love's approximation which omits higher order terras in 
the elasticity relations for the salte of consistency. 

Much work has been done in Soviet Union where 
the mathematical theory of shells is developed by 
Goldenveizer, Novozhilov, Galerkin and others. Vlazov 
has proposed many a simplification suitable for prac- 
tical designs and such theory is generally acknowledged 
as 'engineering theory' of shells. Goldenveizer formulated 
the conditions of computability of strains and showed the 
possibility of identically satisfying the equations of 
equilibrium by means of stress function, which was also 
in similar lines proposed by Lure | 3 |* Vlazov | 4 | expre- 
ssed the internal forces and deformations of a shell by 
two scalar functions. The limit of applicability of 
Vlazov' 3 shallow shell theory are partially elucidated in 
the works of Goldenveizer | 5 |, and Lure | 3 |. 

For translational shells, various numerical inte- 
gration procedures have gained great interest with the 
increasing use of electronic computers. Hedgren and 
Billington | & j combined the stress -displacement rela- 
tions and equilibrium relations into eight first order 
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ordinary linear differential equations. Assuming a levy 
type of solutions for displacements u, v, w, the unknowns 
are solved by Runge-Kutta integration procedure taking 
only one term in the respective series. A good accuracy 
of solution is claimed by taking only a single term. 

Bouma |: 7 | and Apeland |:$;| have proceeded in the same 
lines where they reduce the equillibrium equation to 
eighth order partial differential equation in a stress- 
displacement function. All the unknown stresses and dis- 
placements are given in terms of this function. Mandel 
and Brennan [&\\ used finite difference technique to 
solve for the unknowns u, v, w at each section for a fixed 
edge hyperbolic paraboloid. 

It should be mentioned here that a finite diffe- 
rence technique ia direot but the accuracy of the solution 
depends, to an extent, upon the grid size which is limited 
by the core space of the computer used. Sometimes, to acco- 
modate all finite difference equations, fictitious nodes 
had to be established beyond the shell boundaries. To match 
the number of unknown nodal quantities, some central diffe- 
rence equation had to be replaced by backward difference 
equations. The number of matrix operations to be done are 
also conducive for error propagation. 
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It is also rel event to see to the analysis of 
folded plates by the classical plate theory and elasti- 
city relations. Goldberg and Leve \*Ql developed a solu- 
tion for the stresses in such structures by combining 
the equations of the elassical plate theory and the 
elasticity equations. Applied loading is approximated 
by a Fourier series expansion. Equations are derived 
for the case of a longitudinally simple span structure 
which relate each joint force to a linear combination of 
joint displacements. 4.n number of simultaneous algebraic 
equations are to be solved for each harmonic of Fourier 
expansion where 'n 1 is the number of joints with unknown 
force and displacements. In another paper by Goldberg, 
Glauz and Setlur 1 11> | , the equations-for slabs simply 
suported at their ends - which govern the deformation and 
tractions are written as a system of eight first order 
differential equations. These, together with the joint 
transformation, are integrated using Runge-Kutta fourth 
order process readily giving the displacements, rotations 
and stresses at all points in the antinodal lines. 

This thesis concerns with the analysis of transla' 
tional shell wjhlch is formed by translating an inclined 
line over a parabola. The shell is assumed to be fixed or 
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hinged at the transverse supports and the longitudinal 
boundary may be free, hinged or fixed. The shell is anar 
iy sed considering it to be shallow in the longitudinal 
direction on the assumption that the squares of the 
first order surface derivatives are negligible compared 
with unity. This is also asserted justifyable by E. 
Reissner |12| who states that the shallow shell theory 
will be accurate enough for practical purposes as long as 
the first order surfa.ee derivative is less than or e qual 
to 1/2. 


Basic shell equations are derived in chapter 2. 
The shell equations when shallow in longitudinal direc- 
tion are given in chapter 3, which deals with the funda- 
mental equilibrium equations, in non-dimensionalised 
form. 


Chapter 4 deals with the displacement functions 
corresponding the transverse support boundary conditions. 
Longitudinal boundary conditions are also given in detail 
in this chapter. 

The numerical analysis part of the problem is 
dealt with in chapter 5. The chapter 6 gives the conclu- 
sions drawn from the results. 



CHAPTER 2 


BASIC SHEXL EQUATION'S 

2. 1 Geometry of the Shell ; 

The surface under study is obtained by translat- 
ing an inclined line over a parabola. The equation of 
such a surface in the rectangular cartesian coordinates 
(x, y, z) is given by (see figure 1), 


z 


4 H ^ 

n— ; 

L 


tane.y 


( 2 . 1 , 1 ) 


where H 
0 

L 


height of the shell at x = 0 
inclination of the generating line to the 
y direction 

span of the shell in longitudinal direction 


The position vector to any point is given by, 

R ** xi+yj+zk (2,1,2) 

where I, ic are the unit vectors along x, y, z directions 
respectively. 

a and p are the curvilinear coordinate lines and 
in the sequel subscripts 1 and 2 are used for a and p 
directions respectively. The vectors and eg are ijnit 
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tangent vectors along oc and p directions and e ^ is the 
unit normal vector, 


- e 1* e 2 
e 3 “ AB sinx 


( 2 . 1 . 3 ) 


where 


A = Lame's parameter in ex direction 

B = Lame's parameter in p direction 

x = angle between e ^ and 


This nonorthogonal triad can be eixpressed as, 

z 


'1 


e. 


1 

A 

0 


0 r 

z 


B 




r 

Z 


(2.1.4) 


L 





where z x 


A 

D 


0z „ _ 0z 

Fx » z y Fy 


/ i . 2vl/2 

( ' + a x ) t 


B 


- (1 + Zy)^ 2 - 


/ * , 2 L 2v 1/2 

( 1 + + z y ) 


I 


The egression for the square of the differential 
of arc length, known as the first fundamental forrji of the 
surface, ia given by 
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I = ds 2 ~ A 2 dx 2 + 2 aB cosy dxdy + B 2 dy 2 


z .z 

where Cos X ~ • — t=£ 


The second fundamental form for this surface is given as, 


II ~ hdx 2 -i- 2M clxdy + H dy 2 


(2. 1.5) 


where L " 


M = - -JE = 0, H 


M = o 


also, 


= radius of curvature in a direction =» 


radius of curvature in (3 direction 


= cross -curvature = 


The gauss -weingarten equation are represented by, 


2 - 
9 r: 


( 2 . 1 . 6 ) 


FxFy 


N e. 


Where the value of chri3toffel symbols are given as, 


Wxx 


„ Vxx 

'll _ 

D 


{]! = n ' L . n i . p z - 
• /a I iz ~ I -ix. - / -zjt 
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2.2 Vectors of Elastic Displacement and Elastic Rotation 
oT the Middle Surface : 

Assuming u, v, w are the displacements of a point 
on the middle surface of the shell in a, p and normal 
directions respectively (see figure 2), we can write a 
vector of elastic displacement U of the middle surface of 
the shell such that, 

U = a ^ + v e 2 - w (2,2.1) 

Introducing an auxllliary trihedron such that (see figure 3), 
e^.e| - sin x = ®2*®2 

® 1*®2 “ 0 ’ ® 2*«1 

®1 x ®2 “ sin * 

e. 

i*e, - cot X &2 (2.2,2) 

5 2 “ ill- " 00tX 5 1 

I 

The following notations stand for the strains and rotations 

i 

of the surface, ; 

6 1 = strain in the a direction 

6 2 53 strain in the (3 direction 

<0/ * angle through which the element directed 

along rotates towards eg in the tangent 
plane. 
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= angle through which the element along 
rotates towards in the tangent plane* 

= angle through which the element directed 
along the vector e 1 rotates towards vector 
in the plane (e”, e^). 

^2 “ angle through which the element directed 
along i°wards 1^ in plane (e? 2 » 

co s co, -f- co% = change in angle 'K between the 
coordinate lines. 

6 - ^> 2 - = a measure of the angle of rotation 

of an element of the middle surface around 
from ©2 to e^. 

It can he shown that |30|, 



where « Qtf/ftx* U2 42 9H/9y (2*2,3) 
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We can consider a vector Jz. called tlie vector of elastic 


rotation 

such 


- / 

-a , 

, e 2 

■4. 

• V 

.-a. 

‘ e 3 

i. 9 ® • 

.-/x 




- Y-2. 

6 


= VH ^ + A ^ + 

* SV„ X 


cT^j 


(2.2.4) 


Here the vectors -n- and U are not independent and it can 
be shown that. 


U 


1 


ti » T*" • ® 




^ I 
B • e 3 


u 1 3w 

^ ‘aE 

i aw 
1 ay 


(2.2.5) 


and the components of bending deformation are given by, 


k„ = - t-' ♦ e 


a.; - _ _ . ~x~y a 


* ra.^B, 


7S + ^ fc 


i a / AB\ 

^ D ' 


I 


- 

TT # 1 


i a / ^ AB \ 

B 5y ^ D ' 


. Vy 0 / */ AB\ 

+ SF c ” } 
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‘C - 


-iL t 

A 


5 - - ! L- (^AB) 

9 1 A 3x ' D ' 


T2TT Bx ^ D } * R„ 


A B 


'C i 


•4-z 

"TT 


a. 


25^2,,, a 

„ *JL ® r r, w 

BD By ^ J > 


A 3 0 
15 By” 


(2.2*6) 


where 


0 -a. 


n — y.r^ • 

t B 3T * 


= 


0 -n. 

W 


Introducing new quantities, 




*iAB 



+ 


1 9 
A dx 


( 6 ) 




Jk e 

B * e 3 


1 <L_ 

B dY 


( 6 ) 


2.3 Stress Resultants and Stress Couple Resultants : 

The force vector and moment vector along # - curve 
can he written as (see figure 4* and 5.) 

3?2 ~ Tgi ® “ Tg ®2 ■^23 e 3 (2»3.’1) 

®2 ~ M 2 e 1 " M 21 e 2 
and along P curve, 

T 1 = -T 1 e 1 -'D e 2 + T 15 e ? (2.3.2) 
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f 1 = - M 12 g 1 “ M 1 


where T^ and Tg 35 normal forces/unit length 

*12 , *21 , *13 , *23 = shearing forces/unit length 

and Mg ~ momenta/unit length 
an ^ Mgi ~ twisting moments/unit length 


These forces and moments are expressed in term of normal 
and shearing stresses (see figure 6) as given below, 


T. 


Tr 


+V2 

f , ff; (i + | ) dz 

-h/2 k 2 


+h/2 

f tfl (1 + | ) da 
-h/2 K 1 


*12 ^"uM ^ ^ + f ? ) dB 


+h/2 

f , 

-h/2 

+h/2 

/ , 

-h/2 


*21 " “ /o ^ ^ 1 + S- ^ dz 




-fh/2 


*23 3 -^o + i> d * 


-h/2 


+h/2 

T 15 “ ‘ fh/2 <%l 11 + V ^ 



14 


where 
and ' h 


M. 


+h/2 

/ <r, z (l + |) da 

-h/2 "2 


M, 


+h/2 

/ <r z z( 1 + 4 ) dz 

-h/2 K 1 


M 


12 


+h/2 

/ OTa. z( 1 + ■§ ) dz 
-h/2 R 2 


M 


21 


+h/2 

/ 6- z(l + ~ ) dz 

-h/2 R 1 


( 2 . 5 . 3 ) 


o&i and Otj are normal and ahear stresses respectively 
1 is the thickness of the shell. 

By substituting the values for the stresses we get, 


e, - cotx + 11 e„ 
i = Q ( _J ^ ) 

X ■] ' sin7 / 


sinx 


6o - cotx ^ ^ 6^ 

I, = 0 t ( — ZTTTtc 1 ) 


sin X 


T 12 “ “ T 21 


= o 2 { " - O otx(e 1 + e 2 ) 


sm x 


~ cotx (e.| + 


k. + 21 lc 9 

M 1 = °3 ^ sin'x ^ 
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where 


kn + ■» k, 

Mp = C, ( — ■: ) 

2 3 v sm* / 


*12 3 " 0 


M 


21 


4 sin* 


\ + 


Tr, 


^4 ^ sin / ^ 


(2.3.4) 


~ Poisson’s ratio 
h = thickness of the shell 

„ - Eh 


11 h - f' 2 ) 


°2 “ 


Eh 


2( 1 - v) 


°3 = 


Eh' 


12(1- -u a ) 


°4 = 


Eh 3 

T2n +~v y 


(2, 3.3) 


E = Young’s Modulus of Elasticity 


upon substitution of the values of strain from equations 
(2.2.3 to 2.2.6), we get, 
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T„ = C 


i - 

1 l D 


du 

9x 


w 48 (i + 


T5 


^ C 1 


D 


\ 3u , z x 2 y z xx B 

) u 


A 0v i 
D ByJ 


To = 0 


ii* #V 


(D By 


a x z y z xx 
23^” 


z„z 


B « u - 


AD 


SFo-**^ 


Qx 


A . R 2 2 
AB Z v 2 


R.D 


x z y / „ , V) D N ... , V B 3u l 
rr^ (1 + ^ “2“T ) w + TT 9x a 


z x z y 


T 12= "* 


a , - ",!«!«- 

] 


AB AD B 
2 2 

, 9v f B „ , D L x z y N 

+ 5x KD “ (“2“ * 2" ^ 

ox,i.JJ a^B DAB 


o o 2 2 

Z 2 - A^-r>^ Z Z 

. .. Y XX ( A B L x .y 

+ u t n*~ + if- " u ) 

D IT 


~y 

A^D 


+ w fVl ( B 2 + 2 2) 

a y 


DR 


1 


. z a ( 1 + ^ ) 
5u x y v ' 

" 0x * AD 


ay 

3? 


z z 





17 


M„ = G, 


4 


2 2 
2 . 2. 


3w y xx n , x y N 

57 * z7T ' 1 + ™T! ) 

oy at ir 


z z „ 2 

_xr (1+ « ) . 


nZZ^ n n BzZ ^ 

„ 3 ^ jc xx ^2 + 2A^) - x xx 


R 1 AD 


D 


’Ox 


2 2 2 2 
. B^A 9u , B" Qu , n 3 w 

' H p 57 ' l ” ngr ?rz + -5 “S' 

R^IT 0X IT 0X D 0X 

'B if. 9 2 w _ z x z y A 5u ] 


+ " — S' 

JT 9y 




Kyi 


Mo = G. 


t- D 


2 a 2 w V- A 


3y £ 


R. 


$r 


2 . Qw z y z xx / . . x y ^ 

- ^ 57 “^73 v 1 + ■ 2 ) i 

oy k D ? D 

t. 2 p 

•;j 11 z x z xx /ty 2 , 0 , 2 ^ - 7 , B^ 0 w 

~ R AB 4 " +2A)u ” J Vxx5 * 


2 2 

2,. 2, 


fc r 


, B 2 2J A Du L 7 J B 2 d 2 w 

+ T~ 5x ~2 — 2 

R.j]J 0X D 0w 


i 


M 12= °4 


f - A Vy 8 2 w . a 2 B 9u , AB 8 2 w I 

l 1T7 7^1 
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«21 * 


, , AB 0 2 w B 3 v 3 x 3 XX 0W 

- c 4* 7 5*5y + — J '~4 5y 


B z Z Z„ o o 

+ — VH 2 (3a 2 * 4 ) 

R^D 4 y 

z z B- z z B-*2 
x y 0u x y o w 

‘ T^~ k ' I? - Tx 7 
2 

, z x z y Z xx -d Ow , 1 3u ( 

+ ^5 B • R^Sy) 


2. 4 Equilibrium Equations : 

Let the external loading and external momenta at 
any point be given by, 


P - «i + P 2 ®2 “ *3 


(2.4.1) 


M e = % e 1 + M e2 e 2 


(2.4.2) 


The equation of equilibrium are given by |50|, 

2 

SEi7 lx [ B < T 1 + C03;< I 12 ) j - 1 ( ll/bsin* ® 12 


i a 


sin* 9y 


[a(t 


21 


- cos 


X t 2 )J - B(/J 2 2 ')ainXT 2 


STOT + * Cos * *2> = 


0 
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i a 


sin* S x 


[b(T^ ? + co 3 x T i )] -A ( )ain Ar 


HnT I? t a(t 2 " oosx T 2 J + 5 " ^ ain * 

&r <sf - Sjf> + AB(p 2 + “ a * V = ° 


I. T„ t. - T 


**ij + 4 + ^T 31 ’ + k < B T ia> + fyU * a5 ) 


+ 3I55T Is f B < M i2 + 00£ * M PJ - r 'll sln)! M i 


AB sins P 3=0 9 

. -I ~n £ * 


BIST ly U^ M 2 ■ 


cos M 21 )] *» B ^sinX M 21 


+ AB + AB(M q1 + coax M 02 ) = 0 


-.skr k L b(m i + 003 x M is>] - A 4.' sin * 

+ Sli hr ty [ A(M 21 " ooa * M 2 ~>J- g-fiiSinXM, 


12 

2 


- AB + AB(M 0 2 + cosy M 0l ) = 0, 


p Mp -i Mp “ M 1 

sin 7 (T 12 + T 21 ) + + r~ + _ = o 


(2.4.3) 
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The last of the above equations remains an identity* Upon 
substitution of ohristoffel* s symbols, we obtain, 


8T 1 81 12 2 M 21 A 8T 2 

43 Fx + Vjr TV • A Ty - + B z x z y FV 


.2 D z z 

A- U m + A T) P + X y p =0 
R 1 i 13 AJJ B 2 U 


3 0T 12 ? 3T 1 t? M 21 

4 B TT + z y z xx T 1 + A z y z x FV " r z x z y FV“ 

* 3b *2 + T 2 Vy P 1 " 0 


AP 9T j ^ QTp^r 

+ B “5x + A 7T~- + ^ = 0 




0M, ? 5M* 2 5 Mp a 3M pi 

48 Sx + z x z y 5x A W * 5 Vy Fy 


+ AD T 25 + AD M 0l + *g 2 x z y ^©2 “ ^ 


l3t, «1 


8M 12 . .4 8M 21 


A B FV + z y z xx M 12 + 4 Vy FT" + A 57 

3 a x z v ^2 3 

- A^ -4-2 s— - A^D 


B 


ST “ A D T 13 + DA M e2 


a 2 d 

+ T ¥y M e1 


0 



CHAPTER 3 


SHELL EQUATIONS Y/HEH THE LONGITUDINAL PARABOLIC CURVE 
IS SHALLOW 

3.1 Introduction : 

Most of the shell roof structures come within the 

category of shallow shells where the height to length ratio 

is < l/5. It should he noted that the theory used in a 

shallow shell theory, is based on assumption that the squares 

of the first order derivatives of surface functions, i.e., 

2 2 

z . 2 „z„, are negligible in comparison with unity. This 

y x .y 

limitation on the application of the shallow shell theory 

has been examined by E. Reissner Jiafwho states as a general rule 

that the shallow shell theory will be more than sufficiently 

accurate as long as z , z < l/8 and often accurate enough 

for praotical purposes as long as z , z < l/2. 

jk y — 

2 2 

In the case of shell under study, % , and z have 
been neglected but z z,_ is kept as such since it is not 

x y 

shallow in the transverse direction. The products of l/R 1 
with u and v are neglected in comparison with w in the 
Btrain-displacement relation. Value of sin* is approximated 
to 1 whereas Cos^x is approximated to z„z_ r /B, the derivative 

A y 

of which is z^z^/B which is kept in tact as only the square 
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of surface function derivatives are omitted. Under those 
concessions, we can fairly say that the formulation and 
analysis of the shell under study will be accurate enough 
for values of H/L l/5. 

3.2 Stress Displacement Relations and .Squill ibrium 
Equations ; 

Under the above approximations, the stress dis- 
placement relations (2.3.6) are reduced to, 
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M, = 0, 
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?- -*r 
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xx 9w _ ~x y ^ 
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z v z ( 1+*) 
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0 w 

u xQ y 


Mr 


■\ 2 _ - 2 

. o w , v 0 w 

+ ry +• — tj 1 "n 

0x B* Oy 


afw . V^ 1 + 2l) 

B ta Qy 7 


= °?[7 77 7 


B“ 


~ 2 z z 

, 3 vv ^ 7 xx 

- 1 ' ~ -tar- 07 


5x 


B 


8w £ 

5y ) 


a 2 
0 w 

9x0y 


(3.2.1) 


In the first equillibrium equation, the product 
of b/R.j with is neglected and putting T^g = we 

get, 


dT 1 

B F3T + z - z 


0I 12 Vv 8T 2 

X“y 53T~ + 5jT + m 'i + We “ 0 


0T 12 el 1 Vv 0!C 12 

B 9x * z y z xx + z x z y 9x + B Fy 


+ BP 0 + z x z y P 1 = 0 


•n ^^13 ^23 

I 7 *1 + B W 1 + + ®3 = 0 


9M 12 ,0M 1 QM2 zz 0M 21 

B chic + z x z y S3T “ or + T-^ Sy~ + Bj? 23 


* M e1 * Vy M e2 = 0 
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B 


9M- 9M 1p 0M ?1 

* +ZZ„M ho +ZZ * ^ + * --•-- 

Ox y xx. 12 x y Ox dy 


z x z y ^2 
B Oy" 


BT 13 + M e2 + z x z y M e1 = 0 


( 3 . 2 . 2 ) 


3. 3 Loading on the Structure : 

The loading considered for the structure are, 


a) Uniformly distributed load in plan: 

For uniformly distributed load, the load vector 
is - I? z E where is the magnitude of the load. From 
equation (2.1.4) 


i . e. P ^ 




- P 55 3 

-i- « 


X) 



(3.3.1) 


b) Wind loads: 

The analysis of structure for wind loads is done 
as per 18 Codes \2$\. The loading condition for a single 
fold is given in Table I, For multiple folds, the loading 
remains the same to the windward span and the leeward slope 
of the leeward span as per Table i; On all other roof slopes, 
account need be taken of only the effect of wind drag 
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(see Fig. 7) For wind loads acting normal to the 

longitudinal direction, the load vector can he written 

as, -P T, where P„, is the intensity of wind pressure, 
w w 

From equation (2,1.4), 


P„ = 


_ w 


z y z X A 

~hr- 



P a 
w y 

D 


(3.3.2) 


3. 4 Reduction of Squill ibrium Equations : 

The five equilibrium equations (3.2.2) are 
reduced to three equations in u, v, w directions res- 
pectively by substituting for and T 2 y These equa- 
tions in forces and moments when expressed in displace- 
ments u, v, w, we get, 


-'2 X Sy 2 '**. (-'tg? ^ su .t-e [z~ z ^ £ y 2 Or*? ] A 

^ > Ip. - 2 z,c ^ c ^ % o -■») - 'u V&* ot-»y 

d & d 2Kc2- ' ~jp. 3 ^ 

•* Uh V & -[*’*¥(^-0^1 & - iS »u, 

* v L 6 s - J jeT H 

H- o? ( <- *0 £ . j3_ 


S*/ ^ 


2x_^ 


^ H- *112 - O 


C3. 4-'l) 
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g 2 - 3^2. L g 2 - J dxty 


+ b M-*>; a^r - + f4 - va-7<^pl ^ 

dx 2 - S L B B 3 J 3^, 2 - 

— 2 Z* 2y O y ^ Mr r ^/■^Zy 2 - 71 a ij~ 

o a;* ay l_ ~~^£ J ^ 
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B 3 jp 


g a a a; 3y3 
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The coordinates and displacements are non-dimensional- 


iaed 

as given 

below, 



u = 

u h ; 

v - 

v h ; 

w - w h 

x = 

X L ; 

y = 

r b ; 


tia 

u 

a 1 ; 

u 

a 2 ; 

« a 

ra 

ctf 

11 

v — 

v 4 

where 


Eh'V 12(1 ■ 

X 

D'h 

. v — 


( 5 . 4 * 4 ) 


Using relation (3. 3; 2), the equ.illibriura equations 
are written as, 

*,*,>) 2 + % || + *3 £z + fijr 2jZ 


*x 2 




Jv?- 


+ fl ‘ x 0y + A r * tt * », X 0 +*, * t 

gj? 

2xi y 


My ■ ' '■ + "» " + "■» * syi 


fl„ jv 4 . 

3x2y ^ •— 


fl/iy. X £ - 


o. 


ssr + ^ 

(3-4' 5 J 


E 2 C*, *,*>*.*> 5 6| ^ ^ + A 2x a + * 3yi 


v2- 


t-S.' K^i +6/ lit + S 7 4% 3£ 

^x’ay jx^y 3 a* 5y 

+ e, ^ + Ii„ x 2 £> + b„ x + s„^3a 


+ Bur + B a p + »*. 


3*a ' " a/iy ~"' z ' *y 

8 /Ll x z f = a . 

(3.4.. 6) 
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£3 = 


+ ^ * J§rL + *7 $1. + c 8' * i*£ -I- e<, X aV 




Jx^y 




+ C|6 SyL + C " ^ + C ~ f-t +<=»* ^ 

7 y ^ bx^yl 


+ e /4 x z -f e (b - f - c- 


2 ‘K L dy L 


c*- 4 r 0 


The constants A* s, B’s, C's are given in the Appendix A. 


The stress-displacement relations in non-dimensiona- 


lised form become, 


r 0U 4-n W4.P Y 0 u o.r- 0Y 

G 11 5x + G 12 w + G 1? X QY + G 14 Fy 


G 21 H + G 22 X H + G 23 Ft + G 24 * 


G 31 tl + G 32 §5 + °33 X 5 + G 34 X §S + G 35 X §? 


& + & X SfS + & Ljj + & sfw 

G 41 3¥ + G 42 X OTy . G 43 G 44 ^ 
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V 

D 


r , r y 3 2 W . « 0 2 W r , 9w 

G 51 ^5 + G 52 X 0X9Y + G 53 G 54 By 


M-iotl a 2 a 2** 

12 _ n 9 w , n 8 w 

“5 G 61 + G 62 my 


M 2l b = 

5 


n Q 2 ^ , r 9 2 w 

G 71 OTY + G 72 ^2 


(3.4.8) 


The coefficients G.^ are given in Appendix A. 



CHAPTER 4 


BOUNDARY CONDITIONS 

4.1 Transverse Edge Boundary Co nditions and D is p lacement 
Function s : 

The ’IContorovitch method* is applied for the solu- 
tion of the problem. The functions governing the behaviour, 
in the longitudinal direction, of displacements u, v, w are 
assumed. By the method of virtual work, the equations are 
reduced to ordinary differential equations of unknown 
functions in the transverse direction. Assuming, 

M 

u » 2 u*(Y) f m (X) 

m=1 m m 

M 

* " J 1 4> m W 

M 

w - 2 w*m ^ m (x) 

m- I 

The equilibrium equations are integrated in the 
X-direction by arirtual work principle. 

1/2 

/ E.Cu, v, w, X, Y) 6(u) dX « 0 

- 1/2 1 
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1/2 

f Er,(u, v, w, X, Y) 6(v) cIX = 0 

- 1/2 * 


1/2 

/ E-(u, v, w, X, Y) 6(w) dX = 0 (4*1.1) 

- 1/2 ^ 


Expanding the displacements in terms of displacement 
functions, 


1/2 

f E 1 Cu, v, w, X, Y) 6 U*(Y) f^cx) ax = 0 

t f 


1/2 

/ E 2 (u, V, w, X, Y) 6 v*(Y) dX = 0 

1/ 


1/2 

/ E 5 (u, v, w, X, Y) 6 w#(Y) > B (X) dX = 0 

m = 1 , . . * M (4*1.2) 


Thus giving, 

&ug . E.,(ug, vg, wg, Y) = 0 

6vg . E 2 (ug,. vg, wg, Y) = 0 

6vg , E 3 (ug, vg, wg, Y) - 0 

= 1 


k 


y # • M 


M 


(4*1.3) 
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4 . 2 Stress D is placement Rela t i ons : 


The stress displacement relations are expressed in 
matrix form. 



IP = I? u (4*2.1) 
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4.3 Fixed Boundary Conditi on: 

u = 0 

v = 0 

w = 0 

aw _ n 

n - 0 

The functions assumed are, 

y X ) = 8in(2nmX) m - 1,2, ... M 

</> m C X) = cos(m'jtX) m = 1,3, 0 .. M 

Ym (X) = ( cos(2imxX) - (-1) m ) m = 1,2,... M 

(4.3.1) 

4. 4 Hinged Boundary Condition at Supports : 

At the hinged boundary, 

u = 0 

v = 0 

w = 0 

2 - 

** - — Sr =* 0 
1 ax 

And the displacement functions can be written as, 
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f (X) = sin 2 mu}C 

f m (X) = cosnmx 
y m (X) = cosnutX 


ra * 1 j 2 j a * 4 M 

m - 1 1 3 f ••* H 
m — * 1 , 3 j , h 


(4.4.1) 


Talcing M terms in each of the aeries. 


4.5 Equil i brium Equations and Stress Displacement Relation : 

In the following only one term of the series is taken. 
The accuracy of this has been established in the analysis 
of plates and shallow shells. Taking one term, the equili- 
brium equation become, 


7 * , 7 du* , 7 d_u* - 7 dv» 

V* + A 2 Hy * a 3 * V + A 5 W + 


7 d 2 v* 

a 6 — 2 
0 dY 


_ Hut* m 

+ A^w* + Aq + AgP - 0 


B^u* + 


S + 


5 <i 2 u* 

^ dY 


B^V* + 


“ + g d_vft 


B 


’5 xr 


— ? 
dY 


+ b 7 dw * 


3T + V 


= 0 


® 1 U * + 5 2 w * + 3 3 + 5 4 1? + ®5 $r + 3 6 “rr 


+ 5 o 8 ^4 o 9 p = o 

1 dY^ dY 4 y 


(4*5.1) 
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The coefficients A's, B's, C' s are given in the Appendix A* 
The stress displacement relation become, 

F f = T f u (4.5.2) 

F g = T g u (4.5.3) 

Equations (4.5.2) and (4.5.3) are for fixed and hingfed 
conditions respectively. These are obtained by substitut- 
ing the proper functions f(X), <£(X), Z^(x) in the equa- 

tion (4.2. 1 ). 

4.6 Longitudinal Edge Boundary Conditi on: 

Three main categories of longitudinal edge boundary 
conditions are considered. 

1. End Condition (E.C.): The end condition comprise 
hinged, fixed and free conditions of the longitudi- 
nal edge of the end units. 

2. Symmetric Condition (S.C.): The symmetric conditions 
are formulated so that the integration path can be 
reduced to half provided the structure is symmetric 
with respect "bo geometry, load and end conditions. 

3. Joint Conditions (J.C.)s The joint conditions repre- 

-J-U 

sent the transfer of the variables from an i plate 
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to plate satisfying the comp at ability of deforma- 
tions and forces at the joint. This holds good for 
any ridge or valley of the structure. 

4. 7 End Conditions : 

1) When the longitudinal edges are fixed, 

u = 0 

v = 0 

w = 0 

|f - 0 (4.7.1) 

2) The hinged conditions are, 

u = 0 

v = 0 

w = 0 

m o (4*7.2) 

dY 

3) The free boundary conditions require the forces at 
the boundary to vanish (see figure 8), 

0M ?1 

\ = ^ + S - 0 


0 
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m 2 = 0 

T ^ = inplane shear - T 2 i * M 2 - = 0 


The force boundary condition when expressed in terms of 
displacements are functions of x and y. This has to be 
reduced, to a set of ordinary differential equations 
suitable for numerical integration. 


There are two procedures for doing so. One method 
adopted in to satisfy the end force boundary conditions 
wherever the particular force is maximum. This does not 
ensure the accuracy of solution at other points. But the 
errors introduced by doing so can be kept minimum by 
talcing as many number of terms in the series expansion 
of displacements and satisfying the force boundary condi- 
tions at that many points. 


When only a single or few terras are taken the better 
procedure is to put virtual work boundary conditions* 
Accordintly, 



V n 6(w) dX “ 


0 
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1/2 

/ 1' 6(v) dX = 0 

- 1/2 ^ 


iy z M 2 6( f| ) ax = o 


1/2 

/ T 1 6(u) dX = 0 

- 1/2 1 


Writing these four equations in matrix form, 


B D 1 + E 1) 2 =0 (4.7.3) 


where 



The elements of the matrix E and 'Ef are given in Appendix B* 
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4.8 Symmetric Conditions (S.C. ): 

The symmetric conditions are (see figure 8), 


v^ cos0 i 


w i sin0 i 


0 



0 



= 0 


V = vertical shear - T 0 sin0. + V„ cos0. = 0 

z c. i n i 

If the variation of v and w are not same in the longitu- 
dinal direction, the first equation is satisfied at the 
point where v and w are maximum. The last equation is 
expressed as, 


1/2 

f , v„ 8(w J ax 

- 1/2 


z 


0 


where w = vertical displacement = w. secG. 

3 1 J- 

These four conditions can be written as, 



where 




ft • 


w* . 
m 


m — 1 » o . • M 


( 4 . 8 . 1 ) 
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The elements of the matrices S and S are given in Appendix B. 

4o9 Joint Conditions (J.Q.) : 

The joint conditions at a ridge or valley can be 
expressed as (see figufe 8), 


u i 


U 3 

d 1i 

= 



= 

** 


- 


M 2i 

= 

-rs 

CM 

i 



+ 3 ? v 

** 

= 

p v 

^ 1 i 

=; 

f 13 


The corresponding conditions written in terms of displace 
ments and forces are, 

u i U j 


w ± sin0^ - 

7 i 

cosO^ = 

*3 

cosG^ 

" w 3 

sine^j 

y ± sinD^ + 

"i 

cos8^ = 

S 3 

coa6 . 

J 

+ 7 d 

sinD^j 

, dw N 
( Hr 'i 

- 

/ £w V 




(4.9.1) 



4-1 


I 21 00^ - = - T 23 OOS0J + v n j Sin0 a 

v ni oos6 i + T 2i aln0 i = v nl 00 30 j + T 2j sinS 3 


T 11 “ *1J 


(4.9.2) 


The equation (4.9.2) are modified as, 


£ 1/2 M 2i “ - { 1/2 m 2 3 «t5) i dx 


Hi = T„- 


(4-9.3) 


/ (Tg^ coa ®j[ ** ^ n j_ sin.0^) 6(w^ ainO^ “ v^coa0^)dX 
•— i/2 


= f (~T o .cos0. + V .s.in0 .)6(w. sine. - v„. CO30, )dX 

*^l/2 3 11 J j i l i x 


/ j ( v rL i cos0 i * ^i 9 ^ 11 ^!.) SC^sinQ^ + w^cos0^)dX 


= f (V^cosG . f Tn.sinG.) 6(v.sin0. + w.cos0.)dX 
^ ^ 3 ^3 3 tj J J J 

(4.9.4) 


Equation (4.9.4) can he written as, 
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i.e. 


+ “ 6w U©H 3 > + 

5w !m(©n> + 6v L< e n> = 6 w J m < 0 v 3 > + 


im' w Vj 


®H1 + ©Vi = ©Hj + 0 T3 


©Hi + ®W 


■ ©H3 + ©V3 


(4.9.5) 


Now the equations (4.9. 1 )( 4o 9^ 4 • 9 • 5 ) comprise the joint 
condition. Expressing in matrix form, 


TxJ [D^ = [x 3 ] [D^ 
t A 1il [ D ih + 1^21 ] i “ 

Wjl N 3 + 1^3 J t°2l 3 

The elements of the matrix X^, X^> 
given in the Appendix B. 


(4.9.6) 


A 13’ A 23 are 



CHAPTER 5 


NUMERICAL ANALYSIS 

The numerical integration is carried out essen- 
tially tising the Runge-Kutta Gill procedure, The segmen- 
tation technique is used, where in, the path is devided 
into a suitable number of segments so that the effect of 
a unit change of an unknown variable is felt at the far 
end of the path, thus counteracting the error propagation. 


The functions to be evaluated are devided into 
two groups of variables namely X and T, the correspond- 
ing vectors being and • 


Starting Values for Segment 1: 

Let the starting values be represented by Xc j and 
j (i = 1, . .o n/2; j = 0,1, n/2). In the non-hono- 

geneous solution (j - 0), for fixed hinged and free condi - 
tionp, 



(5.1) 


For the influence coefficient or homogeneous solution 
(3 = 1, ... n/2), each unknown function is given a unit 
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increment. In the case of a free edge where the boundary 
condition is a matrix relation given by equation (4-7-3), 
the unknov/n vector being jpg] , the starting values are 
obtained as, 



I 


■ E ik E K1 hj 


(5.2) 


where I is the unity matrix. 


At the end of the segment 1, the functions are 
divided in two groups namely X.. and T . . (i = 1, ... n/2, 

1 J •** t) 

j = 0,1, ... n/2). influence coefficients are calcula- 
ted ho that the :x group of variables at the end of this 
segment is zero, i.e. 



A 3 + 


0 






(5.5) 


The modified values in T group will he, 


t - 
10 


= T 


io 


+ *13 A d 


(5.4) 


For the homogeneous solutions, the starting values B i . is 
calculated so that the X group variables are unity at 
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the end. of the segment, Thus, 


X ij B 3lc 1 


B- • 
13 


-1 


(5.5) 


The modified T group variables will he, 


= T ik B k3 


13 


(5.6) 


Starting values for Segment 2: 

3?or the non-homogeneous solution, we have, 


io 


= 0 


T 


io 


= ft 1 ) 


'io 


and for homogeneous solution, 



= I 



t 


( 1 ) 

id 


(5.7) 


The starting values for the next segment is 
calculated on the same steps done for the starting values 
of second segment and this procedure is repeated until 
the last segment, except at points where a transformation 
of variables is necessary. 
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Transformation of Variables: 

When a sudden jump in the integration path is 
met with as at the fold lines, the variables are trans- 
formed from the m^* 1 to (m + l)^* 1 segment. The transfor- 
mation relations (4*9.6) can be written as, 


KWl M Mm < 5 - 8 > 

IQ m+1 

-l 

where, [p] = [Xd] H Xy ] 

[Q] = f A 2j) ’[[A-jJ - [A 13 16> ’]] 

[R] = [A^f 1 


The starting values for (m + 



segment are, 



where 


H 

O 

^( ra ) 

io 

X. . = 
ID 

I 

*u = 

ID 

*fc ) - 

d -h (m) 

H ik x ko 

*i" ) = 

Q ik p kj 


.(m) 


( 5 . 9 ) 
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Final Boundary Condition: 

For the hinged or fixed edge boundaries, 


*mo + 


X . q\ 
mo 0 


( 1 ) 


( 5 . 10 ) 


where m denotes the corresponding variable set to aero 

(1 \ 

at the boundary and 1 denotes the influence coeffi- 
cients (correct starting values values of X group) for 
last segment 1. 


For the symmetric condition governed by equation 
(4.Q*1)» we get the corresponding boundary condition as, 

S [*io + hi Oj 1 *] + S Ko + T ij “j 15 ] - 0 

i=1 , . . . m i= 1 , . . n/2 

3=1, . .n/2 

( 5 . 11 ) 

■*lco + G j 1 ^ = 0 k " m, ... n/2, j=1,.,n/2 

And for the free edge boundaries governed by relation 
(4.7.3), 

E [*io + x ii + + hi = 0 

i=1 , . .n/2 
j=1 .n/2 

( 5 . 42 ) 
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From the relations (53.10 to 5.12), the correct starting 
values .for the last segment can he calculated in 

each case. 

Backward Pass: 

In this step, the correct starting values for X and 
1 group of variables are evaluated as, 


G (lc-0 

3 


Die 3 


+ 

J 

lc = 1 , .... 2 


(5.13) 


The relation (5.13) holds good for all segments except 
Where transformation of variables are needed, when, 


r (m) 


r»(m) p 

^jjk ^kl 


n (m+1) 

°1 


(5.14) 


3 


O) s 


is the correct starting values of unknown in the 


first segment. In hinged and fixed edges, the other 
variables are known, where as in free edge, the rest of 
the functions are evaluated by relation (5.2). 


The starting values for the T group of variables 

will be, 


5 (lc+1) 

°3 


t (k) 0 (k+1) (k) 

iD 3 io 


(5.15) 


Having known the starting values of each segment entirely, 
the final integration is done segment wise. 



CHAPTER 6 


NUMERIC A! RESULTS ME CONCLUSIONS 


6 « 1 Intro duct ion 

With this formulation for arched, folded plates, 
where the longitudinal axis is a parabola, some other 
problems can also be considered viz. (a.) parabolic 
cylinder where 0 = 0 (b) conventional folded plate where 

Il/L = 0 (c) plates with small initial curvature. 

This method of analysis is also applicable to the 
above variations under uniformly distributed loads though 
the accuracy of which needs investigation. 

6 . 2 Results and Conclusio ns 

The compilation and execution of the programme for 
the analysis of a single fold with free longitudinal edges 

•j 

takes approximately 6^ minutes in IBM 7044- HO'] system. 
Figures 9 to 23 show the behaviour of a single folded shell 
with fixed transverse supports and figures 24 to 36 repi’e- 
sent shells with hinged supports. 

Effect of H/L ratio 

The value of normal stress T^ depend greatly on 
H/L ratio. The nature of T^ is compressive at all points 
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in the structure for high Il/L ratio. At the central section 
along the ridge, the compressive stress T, increase as Il/L 
ratio is reduced where as at the valley the nature of 1', 
changes from compressive to tensile. At the transverse 
boundary, it is seen that compressive stresses develop at 
the valley and tensile stresses at the ridge portion* The 
normal stress T ^ is found to be higher at the ridge line 
in hinged shells and lower at the valley compared to the 
fixed shells. The bending stresses and Mg increase with 
the decrease in Il/l ratio and this effect is more predomi- 
nant in hinged shells. 

Effect of b/h ratio 

A small increase in is noted along the ridge 
with decrease in b/L, There is a small rise in Tg along 
the ridge. is slightly reduced at the ridge portion 
but remains almost unaffected at the free edge. A reduc- 
tion of Mg is observed with decrease in b/L ratio. 

Ef fect of h/h ratio 

Both in hinged and fixed shells, the bending moments 
M,| and Mg decrease with the decrease in h/L ratio. A rise in 
value of T 1 is observed especially at the ridge portion. The 
value of Tg shows a large Increase in hinged shells compared 


to the fixed case. 
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Effect of e 

The effect of variation of 0 between 36° and 42° 
(which is common for the V type folded plates) is not 
very predominant in normal stresses 1' ^ and T 2 . A small 
increase in and M 0 is noted with increase in 0 in 
both cases. 

Effect of Poisson’s ratio (^ ) 

The variation of u from 0 to 0.15 is found to have 
very little effect on normal stress T^ but the value of 1' 2 
is almost negligible for v = 0 and increases with . The 
value of remains same at the free edge and varies linearly 
at rhe central section with decrease in n . A small increase 
is observed in M 2 with reduction in 0 . 

Discussion of the method 

The function assumed for displacement w in the fixed 
case holds good for Il/L = 0 to h/Ij = 1 /5 * For the hinged 
case, the cosine variation of w may not suit very well the 
geometry of shells with h/l = l/5 to h/L = l/B. For shells 
of low 0 and H/l ratio and high b/li ratio, this method will, 
be most appropriate and is well suited for shallow parabolic 
shells and plates with small initial curvature. 
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